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ASTEROID  DETECTION  AND  MONITORING 
USING  A  SATELLITE-BORNE  VISIBLE  SCANNER 


INTRODUCTION 

Collisions  of  planetary  bodies  with  tbe  Earth  have  played  an  important  role  in  die  geological.  Uological,  and 
climatological  develtqmient  of  the  Earth.  Asteroid  collisions  may  have  had  a  significant  impact  on  biological 
evolution  as  a  result  of  mass  extinctions  brought  on  by  global  climate  diange  [1-3]. 

Over  200  Near  Earth  Asteroids  (NEAs)  have  been  discovered.  The  Aten  and  .^xdlo  class  astermds  have  Earth- 
intersecting  (Mbits.  NEAs  are  usually  smaller  than  S  km  in  diameter  and  have  short  lifetimes  compared  to  the  age  of 
the  solar  system.  They  are  believed  to  arise  fiom  the  decay  of  ctMuetary  nuclei  cm'  fiom  ejection  from  the  main 
asteroid  belt  [4-7]. 

Asteroid  collisions  in  modem  times  have  demonstrated  their  destructive  powo-.  The  1908  Tunguska  meteorite, 
the  1947  Sikhote-Alin  mete(Mite  [8],  and  the  1945  metetM  explosion  over  the  Pacific  Ocean  all  released  energies 
equivalent  to  small  thermtxiuclear  ^vices.  The  recent  near  hit  of  an  asteroid  with  the  Earth  (170,000  km)  [9] 
underscores  die  vulnerability  of  our  planet  to  sudden  devastadoo. 

Based  m  these  concerns.  Congress  ordered  NASA  to  devel(qp  techniques  to  monittM'  for  potential  asteroid 
(Xdlisions.  One  result  of  the  NASA  committee  w(Mk  was  the  identification  of  the  potential  threat  posed  by  asteroids 
as  small  as  SO  m  in  diameter.  However,  objects  that  small  are  not  reliably  detedable  using  ground-based  telescopes 
because  of  atmo^beric  interference  as  well  as  manpower  consideraticMis.  One  of  the  suggesdcms  of  the  NASA 
(XMiunittee  was  to  use  a  passive  detection  system  in  qiace  with  DOD  derived  technologies. 

In  this  report  we  consider  the  feasibility  of  spiue-based,  passive  asteroid  nKMiitoiing  using  Fcxal  Plane  Array 
(FPA)  technology  and  image  processing  tediniques.  In  the  next  secdcm,  we  cmnpare  the  detection  csq[>abilities  cd 
infiai^  and  visiUe  FPAs.  Ba^  on  this  analysis  we  describe  our  proposed  visible  scanner  in  the  third  section. 

ASTEROID  DETECTION 

Asteroid  Detection  in  the  Infrared 

Infrared  imagers  have  Icmg  been  used  to  obtain  higb-resoluticm  data  of  the  Earth's  surface  [10-12].  Image 
processing  algorithms  have  been  develtqied  to  extract  faint  unresolved  moving  targets  fiom  a  tmp^  sequence  of 
infrared  images  [13]. 

The  maturity  of  his  technology  led  us  to  consider  the  possibility  of  using  a  ^pace-based  infrared  (IR)  seamier  to 
monitor  asteroids.  We  assumed  a  detector  noise  etpiivalent  irradianoe  (NEAI)  of  lO*^^  W/cm^  and  a  wavdiand  of  8  to 
14  \un  (tbe  regioa  used  by  Morrison  [14]).  For  an  asteroid  of  radius  R  and  a  distance  r  frcMn  tbe  sensor,  tbe 
irradianoe  at  the  iqietture  in  the  wavebatid  A]  to  A2  is  given  by 

(1) 

/  =  (l-pA)LA(7)dA  , 

Jii 


Manuscript  approved  November  9,  1992. 


1 


SHAFFER,  McHUGH,  AND  WANG 


where  px  is  the  asteroid’s  reflectivity  and  is  the  spectral  blackbody  radiance  at  temperature  T.  Astermd 
temperatures  typically  range  from  180  to  280  K  with  a  spectrally  averaged  temperature  of  255  K  [15].  The 
contribution  firom  die  background  is  negligible. 

With  this  in  mind,  we  solve  Eq.  (1)  for  the  detection  range  r  for  the  cases  I  -  NEAI  and  3  x  NEAI.  The  results  are 
shown  in  Table  1  for  50-m  and  1-km  diameter  asteroids.  For  larger  asteroids,  warning  times  are  proportionately 
larger. 


TaMe  1  —  IR  Detection  of  50-m  and  1-km  Diameter  Asteroids. 


R  (m) 

T  (K) 

r  (km) 

(/  =  NEAI} 

Warning 

(hr) 

r  (km) 

{/  =  3xNEAI) 

Warning 

(hr) 

25 

180 

295 

7.400 

0.7 

4.300 

0.4 

255 

2A74 

21.000 

1.9 

12.000 

1.1 

280 

3.967 

27.000 

2.5 

16.000 

1.5 

500 

180 

295 

150.000 

14. 

86.000 

8.0 

255 

2.414 

420.000 

39. 

240.000 

22. 

280 

3.967 

540.000 

50. 

320.000 

30. 

Although  asteroid  detection  in  the  IR  is  possible,  the  warning  time  (assuming  an  asteroid  radial  velocity  of 
3  km/s)  associated  with  small  asteroids  is  not  sufficient  to  deploy  countermeasures. 

Asteroid  Detection  in  the  Visible 

While  the  emitted  radiation  is  concentrated  at  longwave  IR  wavelengths,  there  is  also  considerable  visible  light 
reflected  from  asteroids.  A  theoretical  ccanparison  of  the  emitted  and  reflected  radiance  of  an  asternd  is  shown  in 
Fig.  1.  The  efficacy  of  using  visible  as  opposed  to  IR  depends  on  the  sensitivities  achievable  in  diese  wavebands. 

The  irradiance  reflected  to  the  sensor  by  a  Sfdierical  asternd  in  the  waveband  A]  to  ^  is  givra  by 

where 

Tg  =  5700  K  is  the  Sun’s  effective  blackbody  temperature, 

Rg  is  the  radius  of  the  Sun, 

fg  is  the  distaiKe  fiom  the  Sun  to  the  asteroid,  and 

Ois  the  angle  between  the  Earth  and  Sun  at  the  asteroid. 

The  NEAI  of  a  sensor  is  defined  as  the  signal  irradiance  that  produces  a  number  of  carriers  Ng  equal  to  tiie  rms 
fluctuation  in  the  number  of  carrios 


N.  =  ANrm. 


(3) 
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X  (^) 

Fig.  1  —  Comparison  of  the  emitted  and  reflected  spectral  radiance  of  an  asteroid.  The  asteroid  is 
assumed  to  be  spherical,  half  lit  (0 s  n/2),  at  temperature  of  280  K  and  reflectivity  of  0.05. 


Hie  number  of  earners  Nj  generated  by  a  signal  flux  incident  on  a  deteettv  of  area  A ^  is  given  by 


N,  =  ,  (4) 

where  Tj  is  sanqile  integration  time  and  t;  is  the  quantum  efficiency  of  the  detector.  Tbe  noise  fluctuations  come 
frcHn  two  indqiendent  sources;  detector  (daik  current)  noise  and  photon  noise.  The  nns  fluctuation  in  die  number  of 
earners  from  detector  noise  is 


AAfi  —  iniuA^TV, 


(5) 


where  is  the  rms  fluctuation  in  the  dark  current  density.  This  follows  the  shot  noise  expression 


irm.  =y2eW 

f 


(6) 


where 

e  is  tbe  electron’s  charge, 

/  is  tbe  total  dark  current  dirough  tbe  detector,  and 
4fis  tbe  bandwidth  of  the  measuremenL 
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The  nns  fluctuatioa  in  die  number  of  carriers  fitxn  {dioton  noise  is 


ANp  *  , 


C7) 


wbere  ^  is  die  badcground  pbolon  flux  on  the  deiecttM’.  To  omnbine  two  indqiendent  randnn  variables,  we  a(M  dieir 
variances: 


anL  =  an}  +an}  . 


(8) 


We  can  now  solve  Eq.  (3)  for  die  noise  equivalent  photon  flux  on  the  detector 


~  ^  +AMdriTI  +  l)  • 


Ibe  NEAI  refened  to  die  senscv’s  aperture  is  given  by 


NEAI  - 


XoAc  X 


(10) 


wbeie 

h  is  Planck's  constant  =  6.63  x  10^  J  s, 
c  is  die  speed  of  light  =  3  x  lO^  m/s^, 

A  is  the  wavelength  of  die  observed  radiatioii. 

is  the  area  of  the  qptics,  and 
Tp  is  the  transmission  of  the  t^cs. 

For  asteroid  detection  using  image  differencing,  we  require  that  the  mmimiim  observaMe  asteroid  radiance  is  3 
times  the  noise  level  or 


I = NEAI 


(11) 


Substimting  this  in  Eq.  (2)  and  solving  for  r  gives 


(cos0+l) 

12 '  'i  1  Arf  he 


j  piLAT,)dX 

it 


(12) 


( 

where  is  given  by  Eq.  (9). 

If  we  assume  the  parameter  values  listed  in  Table  2,  then  Eq.  (12)  reduces  to  r  s  2.9  x  For  a  S0>m 
asteroid,  detection  is  possible  at  a  range  of  7  X  10^  km.  This  cotie^onds  to  a  warning  time  of  28  days  (assuming  a 
relative  asteroid  velocity  of  3  km/s). 
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Table  2  —  Asteroid  Detectioo  Parameteis 


Rf  -  7.0  xlO*  km 

=  1134  cm2 

r,  =  1.5  xlO*  km 

Aj  w  (0.001  in.)2 

o 

II 

»  0.7 

p,  w  0.15 

1?  »  0.9 

X  s  0.75  mm 

*rmr  “  elec/s/(0.001  m.)2 

T(  s  49  ms 

^  =  3.52  xlO^  photon/cm2/s 

PROPOSED  VISIBLE  SCANNER 

We  pnqwse  a  satellite-based  visible  scanning  instniment  that  will  cover  the  entiie  celestial  ^riiere  every  day 
with  a  pixel  size  of  O.QS  miad.  Two  linear  detector  arrays  (bars)  will  produce  a  pair  of  images  with  a  3^-min  time 
delay.  These  two  images  when  processed  will  reveal  asteroids  as  moving  objects  because  cf  parallax. 

The  satellite  will  have  an  altitude  of  770  km  (100  min/orbit).  The  instrument  on  board  the  satellite  will  consist 
of  two  4096  element  linear  arrays  sensitive  in  the  visible  wavelength  region.  Each  detector  will  have  instantaneous 
field  of  view  (IPOV)  of  O.OS  mt^  The  swaths  of  die  arrays  will  be  inclined  by  1 1.49°,  cotrespooding  to  a  3.2-min 
delay  fw  this  orbital  altitude.  The  entire  celestial  S|diere  ctm  be  sanqiled  in  one  day  (IS  orbits).  This  is  dqiicted  in 
Fig.  2. 


SCANS 


Fig.  2  —  Search  pattern  of  the  asteroid  detection  satellite.  The 
detector  array  sweeps  out  a  12*  swath  of  the  celestial  qihere  in  one 
orbit  The  line  of  sight  is  stepped  12*  each  orbit,  so  the  entire 
sphere  is  covered  in  IS  orbits. 


The  delay  between  the  images  from  the  two  bars  wiO  result  in  parallax  motion  of  near  objects  as  illustrated  in 
Fig.  3.  For  objects  1.5  x  10*  km  from  the  satellite,  the  parallax  motion  between  images  will  be  20  pixels 
(assuming  a  O.OS  mrad  IFOV).  Objects  at  3  x  10^  km  distance  will  exhibit  parallax  motion  of  approximately  1 
pixel.  An  asteroid  with  a  dimeter  of  SO  m  will  be  detected  and  tracked  with  this  S3rstem  at  a  distance  up  to 
7  X  10^  km.  Larger  asteroids  will  be  detected  at  correspondingly  larger  distances,  as  shown  in  Fig.  4. 
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Image  from  A 


Fig.  3  —  Use  of  parallax  to  detect  near  Earth  objects,  (a)  An  asteroid  is  seen  at  t  »  0  by 
detector  array  A  and  at  time  T  later  by  array  B.  The  two  arrays  see  it  at  different  anglM 
(depending  on  its  distance)  and  hence  at  a  different  spot  on  the  star  field,  (b)  The  images 
formed  by  the  scans  of  the  two  arrays  show  the  asteroid  in  different  locations  in  the  star 
field.  A  frame  difference  of  the  images  results  in  a  dipole  image  of  the  asteroid. 


I 

F 

E 


Fig.  4  —  Maximum  detection  range  with  a  visible  sensor  as  a  function  of  asteroid  size.  The 
asteroid  is  assumed  fully  lit  with  a  reflectivity  of  O.IS.  The  cutoff  at  3  x  10^  km  occurs 
because  asteroids  more  distant  than  this  will  not  have  enough  parallax  motion  to  survive 
the  proposed  detection  algorithm.  The  dashed  line  indicates  the  detection  range  based  on 
seiuitivity  alone  (Eq.  12). 
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The  dashed  line  is  a  plot  of  Eq.  (12)  for  the  parameters  listed  in  Table  2.  It  represents  the  range  at  which  the 
observed  asteroid  signal  exceeds  3  times  the  noise  level  of  the  detector.  At  this  distance,  asteroids  will  appear  as 
unresolved  point  sources  in  the  scanned  image.  The  solid  line  in  Fig.  4  indicates  the  range  of  detectability  for  an 
asteroid  detected  using  frame  difference  signal  processing.  The  maximum  value  for  the  range  of  detectability  is 
determined  by  the  parallax  angle  between  the  detector  bars.  Larger  asteroids  beyond  this  maximum  range  will  exhibit 
subpixel  parallax  motion  that  will  be  canceled  in  frame  differencing.  This  implies  that  at  the  maximum  parallax 
range  of  3  x  10^  km,  the  observed  signal  from  asteroids  gretder  than  2(X)  m  in  diameter  will  be  much  greater  than 
the  noise  level.  Hiis  could  significandy  improve  change  detection  algorithms. 

Because  of  the  sensitivity  of  a  visible  scanner,  modifications  to  the  present  scheme  would  allow  larger  detection 
ranges  and  longer  warning  times.  The  maximum  parallax  range  in  Bg.  4  could  be  increased  by  increasing  the 
parallax  angle  between  detector  bars.  However,  significant  increases  in  this  angle  could  require  the  use  of  an 
additional  satellite  platfmm.  An  alternate  scheme  could  involve  the  use  of  more  oxnputer  stmage  intensive 
algorithms  that  compare  single  irttages  to  assumed  known  star  charts.  These  modifications,  however,  represent 
improvements  on  an  experimental  concept  oqtable  of  astooid  detection  with  sufficient  warning  time. 

IMAGE  PROCESSING 

NEAs  will  exhibit  parallax  motion  in  the  pair  of  images  obtained  by  this  visible  scanner.  Positive  asteroid 
detection  is  achieved  by  differencing  the  temporal  sequence  of  images.  After  frame  difference  signal  processing 
(FDSP),  fixed  objects  will  subtract  to  zero  (i.e.,  below  the  sensor  noise  level)  in  the  difference  fiame.  Those 
exhibiting  parallax  motion  will  appear  as  dipoles  (positive-negative  pairs). 

The  first  step  in  an  FDSP  is  image  registration  and  resampling.  Image  pairs  are  resampled  on  a  common  grid. 
This  requires  finding  the  displacement  between  them  [16,17]  and  flying  some  sort  of  interpolate  [18,19].  The 
resampled  images  ate  differenced.  Asteoid  detection  is  accomplished  by  identifying  residual  dipoles  in  Ae  difference 
frame.  With  good  registration  and  resampling  algorithms,  residual  clutter  in  Ae  difference  frame  will  be  near  the 
sensor  noise  level.  Therefore,  asteroids  near  the  sensor  noise  level  will  be  detected  as  long  as  they  exhibit  sufficient 
parallax  motion. 

An  important  issue  in  implementing  this  sensor  design  is  the  amount  of  storage  required  for  an  FDSP.  If  dr  is 
the  time  to  collect  a  large  enough  image  to  allow  differencing  and  AT  is  the  time  delay  between  bars  (e.g.,  3.2  min), 
then  the  buffer  size  needs  to  accommodate  Nj(AT  +  2dr)  samples  of  data,  where  Ng  is  the  number  of  samples  per 

minute  collected  by  one  bar.  The  time  increment  At  is  large  enough  to  allow  registration  of  the  partial  images.* 
The  first  Ng  portion  of  the  first  image  is  registered  to  the  second  image. 

The  results  from  the  FDSP  are  downloaded  to  a  ground  station.  Further  processing  on  the  ground  includes 
classification,  orbit  tracking  (using  a  time  sequence  of  difference  frames),  and  Earth  collision  determination  and 
warning. 

CONCLUSIONS 

Congress  has  recently  called  for  the  investigation  of  techniques  to  monitor  fm  potential  asteroid  encounters. 
There  are  unresolved  questions  regarding  the  risk  associated  with  asteroids  of  50  m  to  1  km  in  diameter.  In  the 
previous  analysis,  we  determined  that  a  visible  scanner  would  detect  asteroids  as  small  as  SO  m  in  size  at  ranges  over 
7  X  10^  km.  This  gives  enough  warning  time  for  appropriate  countermeasures.  Infrared  detection,  suggested  by 
other  groups,  is  possible,  but  at  current  sensitivities,  the  detection  ranges  would  be  at  least  400  times  smaller. 

We  propose  a  space-based  visible  scanner  to  monitor  the  entire  sky  every  day.  This  sensor  would  see  asteroids 
in  the  size  range  of  interest  (50  m  to  1  km  and  larger)  at  sufficient  ranges.  Additionally,  the  detection  of  non- 
Earthbound  asteroids  this  small  would  give  new  insight  into  the  distribution  of  asteroid  sizes  and  orbits  and  would 
lead  to  a  better  understanding  of  the  nature  of  the  asteroid  risk. 


*  If  we  assume  a  linear  detector  of  length  4096  pixels,  a  dwell  time  of  50  ms,  an  image  width  of  128  samples,  and 
3  times  oversampling,  the  storage  capacity  (assuming  2  bytes/sample)  required  for  the  two  images  is  100  Mbytes. 
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